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Abstract 

We consider an idealized situation where the Quark-Gluon Plasma (QGP) is described by a perfect, 3+1 dimensional 
fluid dynamic model starting from an initial state and expanding until a final state where freeze-out and/or hadronization 
takes place. We study the entropy production with attention to effects of (i) numerical viscosity, (ii) late stages of flow 
where the Bag Constant and the partonic pressure are becoming similar, (iii) and the consequences of final freeze-out 



04 and constituent quark matter formation. 
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1. Introduction 

Experimental observations of increased collective flow 
in ultra-relativistic heavy ion collisions at CERN SPS and 
BNL RH1C, as well as theoretical considerations [1| sug- 
gest that the Quark-Gluon Plasma (QGP) is nearly a per- 
fect fluid with very small viscosity. Furthermore, it was 
pointed out [2| that shear viscosity is minimal in fluids just 
at the critical point of the phase transition, so measuring 
the minimal viscosity may lead to finding the critical point 
of the phase transition between QGP and hadronic mat- 
ter. Studying the entropy development of QGP systems 
provides an important aid in this. 

A fluid dynamical (FD) description of the nuclear mat- 
ter will be considered here with the example of an Au + 
Au central collision at 65 + 65 A-GeV. Using the MIT Bag 
Model Equation of State (EoS), we consider the physical 
vacuum to have zero energy, except in the volume domain 
where QGP exists. Initially this volume is Vmit- Then the 
matter expands in the flow and reaches a larger freeze-out 
(FO) volume, Vfo- A FD description does not constrain 
the location of the FO. Usually external conditions are 
given to determine it, such as a prescribed FO temper- 
ature or zero-pressure, etc. The fluid dynamical model 
we use [3j, \4\, [5j can therefore run well beyond the FO, so 
that the location of physical FO can be selected afterwards 
as a space-time hyper-surface or layer, where FO occurs. 
In order to select the most realistic FO conditions, the 
temperature and entropy density of the fluid may provide 
important guidance, although the FD model does not use 



these quantities directly, only the pressure and the energy 
density (and eventually the net baryon density) in the EoS. 
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Figure 1: The mean specific entropy, S/N, is shown for three different 
FD computations (N is the number of participants). Although the 
simulations were done for adiabatic expansion of an ideal fluid, the 
entropy increases due to the numerical viscosity of the method. The 
difference in initial specific entropy between the two cases describing 
collisions with impact parameter b = is due to coarse graining. 
V ce ll denotes the cell size of the computational grid. 

The 3+1 dimensional, FD model [3j, |4J, [5j uses the Par- 
ticle in Cell (PIC) method adapted to ultra-relativistic 
heavy ion collisions. In his method, so called marker parti- 
cles, which are smaller Lagrangian fluid cells corresponding 
to fixed baryon charge, move in a grid of larger Eulerian 
cells. The computational fluid dynamical (CFD) calcula- 
tion, describing the participant zone of the reaction con- 
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taining a net baryon charge N, starts from an analytic 
initial state model jg, |9( based on longitudinally expand- 
ing strings of the color-magnetic field allowing for different 
string tensions, string cross-sections and for finite impact 
parameters. The entropy of the initial state depends on 
the resolution of the spatial grid, as the coarse graining 
with the relatively small initial spatial size leads to an ob- 
servable difference (see figure Q]). The initial entropy is 
substantial, it is in agreement with other estimates, e.g. 
[6j. The subsequent expansion shows a weak entropy in- 
crease, 5-6 %, due to the numerical viscosity, although 
the model considers a perfect fluid. The entropy increase 
due to numerical viscosity is smaller when the cell size is 
smaller. At late stages the entropy increase is weaker due 
to the smaller gradients between neighbouring cells in the 
expanding volume. The entropy increase due to numer- 
ical viscosity exceeds the one obtained in 1+1 and 2+1 
dimensional models with smaller cell sizes [7[, and it is 
somewhat below the entropy increase of 6-24 %, obtained 
in those models due to physical viscosity. 

The MIT Bag Model EoS yields a negative pressure 
at low energy densities (which occur at late stages of the 
expansion) indicating a tendency for clusterization. As in 
most CFD approaches, the pressure is cut at zero, i.e. no 
negative pressure is allowed in the numerical calculation. 
This is necessary for the stable calculation of the late ex- 
pansion stages. This choice also ensures that the regions 
past the FO surface, most of which will have zero pressure 
in the FD calculation, have a minimal influence on the re- 
gions behind it, as these cells exercise no work, pdV , on 
their neighbours. 

The MIT Bag Model assumes two components: an 
ideal parton gas and a perturbative vacuum. The pertur- 
bative vacuum has a vacuum energy, BV , where B is the 
bag constant (in our calculations B — 0.396 GeV/fm 3 ), 
and V is the volume of the QGP. Thus, this energy is 
increasing as the plasma expands. 

This energy increase happens at the expense of the 
work done by the partonic pressure, p p dV, which partly 
accelerates the neighbouring cells (causing the observed 
flow), and is partly invested into the increased energy of 
the Bag. It is assumed that the parton gas is an ideal 
Stefan-Boltzmann gas, where p p = e p /3 so the Bag model 
EoS with e = e p + B and p — p p — B becomes: p = 
e/3 - 4B/3. 

2. Final stages of the fluid dynamical expansion 

Since the late stages of the expansion are important for 
the subsequent description of hadronization and freeze-out, 
we scrutinize the temperature and entropy development at 
these stages in FD. For transparency, we use the MIT Bag 
Model EoS, and compare it qualitatively to more realistic 
EoS features from Lattice QCD. 

When the total pressure drops to zero then the total 
pdV is also zero. This usually happens in surface cells, 
which continue to participate in the expansion due to their 



inertia. At the same time, if the energy density of the 
Bag remains constant, B, the total energy of the Bag in- 
creases further. It must be noted that all these calcula- 
tions were done with the assumption that the Bag field 
is uniform and homogeneous, and thus its entropy van- 
ishes. Obviously, near the hadronization and FO the ap- 
plicability of the MIT Bag Model EoS is overstreched. The 
Bag field should gradually disappear and hadrons will be 
formed. The assumption of an equilibrium phase mixture 



111 . |12[ leads to a slow hadronization, which is not sup- 



ported by the spacetime extent indicated by two particle 
correlations. An alternative supercooling followed by fast 
and simultaneous hadronization and FO is advocated for 
in refs. [LJ, |l4j, [l5j. We apply this approach so that the 
FD stage of the model concludes in the final simultaneous 
hadronization and freeze-out. 

We shall study the details of this process. To have a 
better insight we shall consider the Bag (i.e. background 
field energy) and the parton gas as two separate compo- 
nents of a coupled system. 

To derive the thermodynamic properties of the QGP, 
such as temperature or entropy density, it is necessary to 
decompose the total energy density e provided by the FD 
calculation into the energy density of the parton gas (e p ) 
and the energy density of the Bag (es): e = e p + es. Then 
thermodynamic quantities can be obtained from the EoS 
of the parton gas. Doing this decomposition is trivial when 
the pressure if strictly positive and es — B is a constant. 
However, when p = 0, several choices are possible: Since 
the system is described as consisting of two parts, and can- 
not be treated as a single thermodynamic system in equi- 
librium, the decomposition can only be done unambigu- 
ously if (1) a particular change in state is assumed and (2) 
explicit assumptions are made about the non-mechanical 
energy exchange between the two components (the Bag 
field and the parton gas). 

Since the FD calculation assumes an ideal fluid, the 
change in state is ideally adiabatic for the whole (coupled) 
system, i.e. once the pressure p has reached zero, the work 
done on one Lagrangian cell by its neighbours vanishes, 
pdV = 0, and the energy E = E p + Eb of a Lagrangian 
fluid cell will stay constant during the expansion. For the 
energy exchange between the two sub-systems within one 
cell (Bag field and parton gas) several choices are possible 
when the pressure becomes zero. For simplicity, let us 
study a baryon-free Stefan-Boltzmann gas, and consider 
the following possibilities: 

The first possibility is assuming that the energy den- 
sity, &B-, of the Bag field stays constant during the ex- 
pansion, es = B = const. The energy of the Bag field, 
BV, increases strongly during the expansion. This is only 
possible if energy is transferred non-mechanically (i.e. in 
another way than through the work done by the partonic 
pressure p p ) from the parton gas component to the Bag, 
and thus the entropy of the parton gas decreases. 

A second possibility is assuming that the only energy 
exchange between the two sub-systems is due to mechani- 
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Figure 2: The energy density of the Bag field, eg, as a function of 
the temperature of the partonic gas, T, for adiabatic expansion (full 
line, dS = 0), constant energy of the Bag field (dashed line), and 
without change of the energy of the parton gas (dotted line). Case 
dE p = assumes that the total energy of the ideal gas component 
remains constant, thus the total energy of the Bag field remains 
constant as well during the final expansion of the system, while its 
energy density decreases. Realistically, the Bag field must disappear 
in the transition, so the real trajectory must lay below the dotted 
line. This disappearance of the Bag field may start before the total 
pressure drops to zero. 



cal interaction (the work p p dV, and thus the entropy of 
the parton gas is constant, dS p = 0. In an adiabatic 
expansion of the Stefan-Boltzmann gas starting from Vo, 
the energy density decreases e p /e p o = (V/Vo)~ 4 ^ 3 : where 
epo = 3.B when the total pressure reaches zero. As the 
total energy is conserved {e p + es)V — (e p0 + B)Vq, the 
volume ratio can be expressed via the energy density ratio, 
so (e p {T)/3B) = [(e p (T) + e B (T))/(4B)] 4/3 . From here 
eB (T) can be expressed as shown in figure [2] 

The third choice considered here is that the energy of 
the parton gas in a Lagrangian cell stays constant after 
the pressure has reached zero: E p = const., or dE p = 0. 
This is possible if in addition to the mechanical work that 
the parton gas does on the Bag field, it also draws energy 
from it non-mechanically, resulting in a zero net energy 
exchange. Thus the entropy of the parton gas is increasing, 
i.e. this is a dissipative expansion. 

Finally, if the net energy transfer is directed from the 
Bag field to the parton gas, then the entropy increase is 
even larger. This possibility is the physically most realistic 
one. 

For these cases, the energy density, es, of the Bag field 
is plotted as a function of the temperature of the parton 
gas for all three cases in figure [5] 

For the purpose of comparison with Lattice QCD sim- 
ulations, the interaction measure, im = (?ic) 3 (e — 3p)/T 4 , 
can also be computed for these three cases. The results 
are plotted in figure |31 and compared with the prediction 
of a lattice QCD calculation. [101 ] 

In the numerical FD model after about 3-4 fm/c time 
some fluid cells reach the zero pressure, and as time pro- 



Figure 3: The interaction measure, IM, as a function of temperature 
for the same three cases as in fig. [2] These curves are compared with 
the interaction measure obtained from a lattice QCD calculation, 
represented by empty circles. There is relatively good agreement 
between the prediction of the Bag Model and that of the lattice 
QCD calculations above a temperature of 200 MeV. [1(]| In a more 
realistic situation when the Bag field diappears in the final expansion, 
the trajectory would lay below the dotted line. 



gresses the number of such cells is increasing. The zero 
pressure cells appear early in the calculation because the 
initial state assumes a smooth surface with low density for 
cells at the surface of the fluid. Although the FD model 
development is only governed by the pressure, the entropy 
and the temperature development after this time depends 
on the assumptions of the energy transfer between the ho- 
mogeneous Bag field and the parton gas. According to the 
first choice the entropy of the system may remain constant 
and then energy should be transferred from the Bag field 
to the parton gas. 

A FD simulation was run for the collision of two Au 
nuclei at 65 + 65 A-GeV energy and impact parameter 
b = 0. The entropy evolution of the system is plotted for 
all three cases (see figure H]). The entropy increase in the 
first case (solid line) is due to the numerical viscosity of 
the CFD model. 

The es = const, choice leads to an entropy decrease. 
The second law of thermodynamics rules out the possibility 
of the entropy decreasing, therefore the energy density of 
the Bag field cannot be constant. More generally, energy 
cannot be transferred non-mechanically from the parton 
gas to the unifrom and homogeneous Bag field, and so all 
such possibilities with total entropy decrease are excluded. 
However, energy transfer in the other direction is allowed, 
as this corresponds to an increase in entropy. 

The final hadronization and freeze-out may well pro- 
ceed through a Quarkyonic matter phase [16J, [l7[ , which 
suggests that first the deconfinement ceases (i.e. the Bag 
field fragments in our picture) , and then the chiral symme- 
try will be broken, so that the current quarks gain mass 
and become constituent quarks. It is suggested [3j| that 
the Constituent Quark Scaling 18|, ll9| may provide a pos- 



sibility to extract information from the Quarkyonic Matter 
state. 
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Figure 4: Results for an Au + Au collision at 65 + 65 A-GeV energy 
at impact parameter b = 0, from a CFD calculation with the Particle 
in Cell (PiC) method with cell size dx = dy = dz = 0.575 fm. The 
mean specific entropy of the Au + Au system, S/N, as a function 
of time in the numerical fluid dynamics simulation of a heavy ion 
collision. Solid line: adiabatic expansion of the ideal gas component, 
dashed line: eg = B = const., dotted line: E v = const. The slight 
entropy increase in the "adiabatic" case is due to numerical viscosity. 

The present calculation indicates that in the dynamics 
of this development and the disappearance of the Bag field 
the entropy plays an important role, especially the entropy 
arising from the structure of the breaking up Bag field. 

Due to the interaction of gluons, we can anticipate that 
in a given domain of space, the energy density of the Bag 
may only be e_e = B (if the field is present), or e# = (if it 
is not). Thus an average decrease in the energy density of 
the Bag field, es, corresponds to the field breaking up into 
smaller pieces. When this happens, the field will not be 
uniform and homogeneous any more, therefore its entropy 
will need to be taken into account. The present study 
provides valuable guidance, indicating the possibilities of 
the detailed description of the break-up of the background 
field and the subsequent breaking of the chiral symmetry. 
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